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We present experimental results on the two-body loss rates 
in a magneto-optical trap of metastable helium atoms. Abso- 
lute rates are measured in a systematic way for several laser 
detunings ranging from -5 to -30 MHz and at different inten- 
sities, by monitoring the decay of the trap fluorescence. The 
dependence of the two-body loss rate coefficient /3 on the ex- 
cited state (2^P2) and metastable state (2"^ Si) populations is 
also investigated. From these results we infer a rather uniform 
rate constant Ksp — (f±0.4) x 10~^ cm^/s. 

PACS 32.80. PjOptical coohng of atoms;trapping, 34.50.Rk- 
Laser modified scattering and reactions 



I. INTRODUCTION 

Helium atoms in the metastable triplet state 2^ S\ (He*) 
appear to be a good candidate for Bose-Einstein Condensa- 
tion (BEC) according to theoretical predictions Q . The cross 
section for elastic collisions between spin-polarised metastable 
helium atoms is expected to be large, allowing efficient ther- 
malization and evaporation in a magnetostatic trap, which 
is the standard technique to reach BEC [^[j^. On the other 
hand, very high autoionization rates (Penning collisions) pre- 
vent reaching high densities of metastable helium atoms, both 
in the presence and in the absence of light, unless the sample 
is spin polarized. 

If a metastable helium atom collide either with an other 
metastable atom, or with an helium atom excited in the 2^P2 
state, the quasi molecule formed can autoionize according to 
the following reactions: 



He (2^51) + He (2^ Si) 



He(2^P2)+He (2=^ Si) 



He(l^S'o) + He+ + e" 



He+ -I- e" 



He(1^5'o) + He+ + e" 



HeJ -I- e" 



(1) 



(2) 



A first experiment at subthermal energy (E = f.6 meV) 
with the metastable helium system was performed by Miiller 
et al. Q , allowing the determination of the interaction poten- 
tials. Using those potentials the rate f5ss for the reactions 
(0) has been calculated to be a few 10"^" cm^/s, which 
agrees with measurements performed in Magneto-Optical 
Traps (MOT) (io[0jl^ • According to theoretical predictions 
[0, the ionization rate corresponding to the reactions (^) 
should be suppressed by four orders of magnitude in a magne- 
tostatic trap. Spin polarization of the atoms and spin conser- 
vation in the coUisional process are the causes of this suppres- 
sion, which makes the quest of BEC reasonable. Actually, a 
reduction of more than a factor of 20 in the two-body loss rate 



in an optically polarized sample was observed experimentally 

In presence of light exciting the transition 2'^ Si 2^P2 , the 
reaction ^ is dominant. "Optical collisions" with metastable 
helium atoms were measured to have surprisingly large cross 
sections when compared with alkali systems The study 

of optical collisions is of fundamental importance in order to 
optimize the first step towards BEC, consisting in pre-cooling 
and trapping the atoms in a MOT. The goal is to transfer 
a cloud as dense as possible in a magnetic trap, in order to 
increase the elastic collision rate and start evaporation. The 
experimental study of optical collisions is the subject of this 
paper. 

Several groups reported measurements of optical collisions 
rates, by studying losses in the MOT at small detunings 
[|o[||0 around -5 MHz and at large detunings at -35 
MHz and —45 MHz. Measurements over a broad range of 
detunings, from —5 MHz to —20 MHz, were reported in 
and the dependence of the loss rate on the intensity of the 
MOT laser beams was investigated. In reference Q a theo- 
retical model for optical collisions is also proposed predicting 
rates in good agreement with the measurements, but differ- 
ing by more than one order of magnitude with all the other 
measurements previously quoted. 

Our measurements are performed in a MOT loaded with 
f 0^ atoms, at a peak density of 10^*^ atoms/cm'^ . With respect 
to previous works, we extend the measurements of the two- 
body loss rate to a wider range of detunings and intensities 
with a good precision, by measuring the number of atoms 
and the size of the trap using absorption techniques. Also, 
by measuring accurately the excited state population in each 
trapping condition, we are able to interpret our data with a 
simple model, expressing the two-body loss rate in terms of 
the excited state population and of a rate constant Ksp, found 
to be independent of the laser detuning and intensity. 

Our experimental setup is described in section Q while in 
section 



HI 



we explain our detection system and we give the 
working conditions and performance of our magneto-optical 
trap. In section ^ we describe in detail the experimental 
procedure used to measure the two-body loss rate and the ex- 
cited state population for different trapping conditions. The 
results are given in section and the conclusions in section 

0- 



II. EXPERIMENTAL SET-UP 

A beam of metastable helium atoms is generated by a con- 
tinuous high voltage discharge in helium gas, cooled to liquid 
nitrogen temperature. Radiation pressure on the metastable 
beam allows one to increase its brightness, and to deflect it 
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from the ground state hehum beam The metastable 

atoms are then decelerated by the Zeeman slowing technique 
and loaded in a magneto-optical trap (MOT) in a quartz cell 
at a background pressure of 5 x 10~^° torr. More details on the 
experimental setup will be given in a forthcoming paper |^ . 
MOT parameters for optimal loading of the trap are listed 
in table |. For the laser manipulation of the atoms, we use 
the line at 1083 nm, connecting the metastable triplet state 
2^5*1 to the radiative state 2^P2. The saturation intensity hat 
for this transition is 0,16 mW/cm'^ and the linewidth r/27r 
is 1.6 MHz. Our laser system consists of a DBR laser diode 
(SDL-6702-H1) in an extended cavity configuration, injecting 
a commercial Ytterbium doped fiber amplifier (IRE-POLUS 
Group). The diode is stabilized by saturation spectroscopy 
at -240 MHz from resonance. At the fiber output we obtain 
600 mW of power, in a TEMOO mode at the same frequency. 
The estimated linewidth is around 300 kHz. All the frequen- 
cies required for coUimation, deflection, trapping and probing 
are generated by acousto-optical modulators in a double pass 
configuration, while we use directly part of the fiber output 
beam for slowing the atoms. 

TABLE I. Optimal loading parameters of the He* mag- 
neto-optical trap. 



Laser detuning 


-45MHz 


Laser beam diameter 


2 cm 


Vertical laser intensity (Ox) 


2x9 mW/cm^ 


Longitudinal laser intensity (Oy) 


2x9 mW/cm-' 


Transverse laser intensity (Oz) 


2x7 mW/cm^ 


Total intensity 


50 mW/cm^ 


Weak axis magnetic field gradient 


hx = by = 20 G/cm 


Strong axis magnetic field gradient 


= 40 G/cm 



III. DETECTION SYSTEM AND 
CHARACTERIZATION OF THE MOT 



Fluorescence PD2 



Absorption PDl 




CCD camera 
for absorption 



A. 



FIG. 1. Detection set-up. By rotating the A/2 plate, one can 
create either a progressive plane wave for measuring the absorption 
on the photodiode PDl, or a standing wave, with both beams cir- 
cularly polarized in the cell region, for imaging the cloud onto the 
CCD camera. PD2 monitors the fluorescence of the MOT. 



In order to fully characterize the cloud, we use a probe laser 
beam on resonance, whose diameter is about 1 cm, which is 
turned on 100 ^s after the MOT field and light beams have 
been turned ofi^. Our detection setup (see fig. |^) allows dif- 
ferent measurements. With the combination of A/2 plates 
and polarization beam splitter cubes, we can create either (i) 
a progressive wave, circularly polarized, passing through the 
atomic cloud towards a photodiode (PDl in figure |^), giving 
the total absorption by the atoms, or (ii) a stationary wave, 
also circularly polarized, one arm of which is sent to a CCD 
camera, allowing spatially resolved absorption pictures of the 
cloud. A second photodiode (PD2 in figure 0) is used to col- 
lect the cloud fluorescence. We use the absorption photodiode 
PDl to measure N, the number of atoms in the steady state of 
the MOT. The probe beam saturates the transition when the 
incident power exceeds 10 mW (see fig. The maximum 
absorbed power is then P — Nhv^. Our Watt-meter (Co- 
herent lab-master) is calibrated to 3% accuracy and allows a 
rectilinear calibration fit of the photodiode voltage. We mea- 
sure a maximum total absorption of 1 mW, corresponding to 
(IzbO.l) X lO" atoms. We estimate the accuracy for the mea- 
surement of N to be about 10%. 
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FIG. 2. Absorbed power by the MOT versus incident power of 
the laser probe beam. The absorbed power saturates at 1 mW for 
an incident power of 10 mW. The corresponding number of atoms 
is (litO.l) X 10'' atoms 

The typical parameters of our magneto-optical trag with 
the operating conditions of table | are listed in table |" 

TABLE II. Characterization of the MOT with parameters 
of table Q. 



Number of atoms 


TV = (liO.l) X 10"* 


RMS size (weak axis) 


(Jx = (yy ~ (2±0.1) mm 


RMS size (strong axis) 


= (1.6±0.1) mm 


Density at the center 


(1±0.25) X 10"' atoms/cm^ 


Temperature 


1 niK 



2 



We stress the fact that the case of He* differs of that of 
alkahs, for which the imaging method gives a direct measure- 
ment of both the two-dimensional column density and the rms 
sizes of the MOT, by absorption of a brief and low intensity 
probe pulse (/ <^ Isat)- In the case of He*, the quantum 
efficiency of the CCD camera (10~^ at 1.083 fj.m) is too low 
to provide images with a sufScient signal to noise ratio. We 
need instead to illuminate the atoms with a 200 fis pulse 
whose intensity is about 0.1 mW/cm (/ ~ Isat), and use a 
moderate magnification of 1/5. Another difficulty with He* 
occurs from the large recoil momentum fik/ra (9.2 cm/s) due 
to the light mass of the atoms : the atoms are pushed out of 
resonance during the 200 /xs pulse if a traveling wave pulse is 
used. The solution we adopted is to illuminate the atoms in a 
standing wave with the set-up shown in figure ^ , Though this 
scheme allows us to obtain pictures with a good contrast, the 
drawback is that the images obtained in the standing wave 
configuration for / ~ Isat are more difficult to analyze than 
in the low intensity case. In order to interpret the absorption 
pictures in the standing wave configuration, and for any sat- 
uration parameter, we developed a handy theoretical model 
(see appendix ^ giving the column density of the atoms for 
each pixel of the CCD camera. The resulting density is then 
fitted by a Gaussian curve to extract the size of the cloud. 



IV. MEASURING THE TRAP DECAY BY 
FLUORESCENCE 

Once the loading of the MOT is interrupted, the evolution 
of the number of trapped atoms A'^ is given by the following 
equation: 



^ = -aiV-/3 / n\v,t)d 



(3) 



where n(r, t) is the atomic density at position r and time t, 
a is the decay rate due to collisions between trapped atoms 
and the residual gas, and 13 is the two body intra-MOT loss 
factor. Assuming that the spatial distribution is independent 
of the time evolution of the number of atoms, which is valid 
at low enough densities, one can write the density as 



n(r,t) 



N{t) 



{2tv)^ axUyaz 



(4) 



At low enough pressure and high enough density, losses due 
to background gas are negligible, so that the equation reduces 
to 



dN 
H 

whose solution is 

N{t) = 



(47r) 2 (7a:(Jyaz 

N{to) 



1 + ^n{0,to){t-to) 



(5) 



(6) 



where to is the initial time. In order to follow the evolution of 
the number of trapped atoms, we monitor the fluorescence de- 
cay of the MOT with a photodiode (PD2 in figure |l]). As the 



fluorescence signal is proportional to the number of atoms, we 
obtain a fluorescence decay curve reproducing equation 
which we fit to get the parameter /3n(0,to)- In order to de- 
termine P, one still has to measure n(0, to), which means that 
one has to measure the rms size of the cloud along the three 
directions and the initial number of atoms N{to). 
Our goal is to measure the loss rate for a wide range of de- 
tunings and intensities. The experimental procedure, divided 
in three successive steps, is the following. 
(1) First, we load the trap for 1 s at 5 = —45 MHz and at the 
highest intensity in the trapping beams (I/Isat=50 per laser 
arm). Then, we stop the loading by blocking the slowing 
beam with a mechanical shutter. 20 ms later, we "compress" 
the MOT by suddenly changing its detuning and intensity 
using acousto-optical modulators. We record the fluorescence 
signal during this procedure. A typical fluorescence curve is 
shown in figure ^. 
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FIG. 3. Evolution of the fluorescence signal. Once the loading 
is stopped, scattered light from the slowing beam is blocked, which 
explains the drop of the signal at t=-10 ms. The detuning is then 
set to S = —20 MHz at t=0 ms and the fluorescence decays. 

The loading is stopped at t=-20 ms and the photodiode sig- 
nal drops by a factor of 2 at t=-10 ms because the background 
light from the slowing beam is blocked. The fiuorescence is 
greatly enhanced in the beginning of the compression phase 
at t=0 ms, as expected when the detuning is set closer to res- 
onance (the detuning is set here to -20 MHz), but decays to 
almost zero in about 100 ms because of the two-body losses. 
Figure ^ shows the time evolution of the size of the cloud dur- 
ing this phase of compression, showing that 10 ms are enough 
to reach the new equilibrium size. Thus, we extract the pa- 
rameter /3n(0,to) from a fit of the fluorescence decay starting 
from t — to = 10 ms. At this very time we measure the sizes 
of the MOT along x and y and the number of atoms in order 
to calculate n{0,to). 
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FIG. 4. Size of the MOT during the compression phase. The 
new equilibrium is reached after 10 ms 

(2) Then, the sizes along the weak axis of the magnetic field 
gradient are measured by absorption on the CCD camera as 
explained in section [II , Figure H shows the rms size along x 
for various laser detunings and intensities. The size along z 
(strong axis of the quadrupole field) is inferred from measure- 
ments of the sizes along x and y with a magnetic field gradient 
b twice as large. We find a typical size along z 20% smaller 
than along the weak axes of the quadrupole. We did not cor- 
rect the sizes for the expansion of the cloud during the pulse 
lasting 200 ^a, as this would have required the measurement 
of the temperature for all the detunings and intensities. Nev- 
ertheless, we performed some time of fiight measurements, 
giving temperatures ranging from 0.3 mK at -10 MHz to 1 
mK at -40 MHz, from which we estimate that the sizes are 
overestimated at most by 5 % at -25 MHz and by 15 % at -5 
MHz. In addition, we measured the statistical error on the 
sizes to be relatively small at large detunings, 2 to 3%, but 
larger at small detunings (about 10% at -5 MHz). This is due 
(i) to the poor spatial resolution of our imaging system (pixel 
dimension x 130/im), and (ii) to a low signal to noise 

ratio for small detunings where the loss rate is larger, as most 
of the atoms are lost during the compression phase. 
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FIG. 5. Rms size of the MOT cloud as a function of the intensity 
of the MOT laser beams for various detunings. 



(3) Finally, to determine the number of atoms that were 
still trapped at to = 10 ms, we simultaneously switch off the 
magnetic field and set the trapping beams on resonance at 
to, instead of letting the trap decay as in figure ^. The laser 
intensity is set to a high enough value to strongly saturate the 
transition. We get a peak of fiuorescence, whose amplitude is 
proportional to the number of atoms. We compare it with the 
peak obtained with the same procedure but for the MOT in 
the best loading conditions of figure ^, for which we measured 
the number of atoms precisely. From this comparison, we infer 
the number of atoms at t = to in the compressed MOT, and 
thus determine n(0, to). 

This measurement also gives access to the value of the average 
population of the excited state Tip. Indeed, -Up is given by 



F 



F„. 



"P n 

Y7^ = 2x., 



(7) 



where F is the fluorescence signal we measure in the com- 
pressed MOT at to, and Fj^^ax the fiuorescence signal at res- 
onance, when the transition is saturated, and vTp expected to 
be 1/2. 

Figure ^ shows the results of the fiuorescence measurements, 
giving Fmax/F as a function of the inverse of the laser inten- 
sity I for various detunings. It is interesting to note that the 
inverse of F is found to vary linearly with the inverse of 7. 
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FIG. 6. Fluorescence signal F from the MOT as a function of 
intensity I of the laser beams. The inverse of the fluorescence F 
is found to vary linearly with the inverse of the intensity /. The 
results are used for the calibration of the number of atoms. 

Following jl^, the fiuorescence of N atoms in the com- 
pressed MOT can be modeled by the following equation: 



F ^TjNhu- 



l-fCaj^^-f 4 



(8) 



where rj is the detection efficiency, I is the total intensity of 
the six MOT beams, and Ci and C2 phenomenological factors. 
Ci and C2 would be 1 for a two-level atom, but they are 
expected to be smaller for an atom placed at the intersection 
of 6 differently polarized laser beams, as happens in a MOT. 
In reference m/A , Ci and C2 are found to be equal, and slighly 
larger than the average of the squares of the Clebsch-Gordan 
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coefficients over all possible transitions. For a, J — 1 < > J = 

2 transition, this average is 0.56. We can rewrite equation (| 
as 



F ^Ch Ci F 



(9) 



where -Fmaa: = r)Nhv^. 

The results of figure g show a good agreement with (^). But, 
C2 and Ci are not found equal, and both depend on the de- 
tuning. For example, Ci is found to be 0.58, 0.48, 0.46, 0.44, 
0.22 for 5 =-25, -20, -15, -10, -5 MHz respectively. We stress 
the fact that, for the fluorescence at resonance, and for full 
saturation, CI and C2 are expected to be equal. 



V. RESULTS 

The results of the Penning collisions rate /3 are shown in 
figures 1^ and ^. 
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FIG. 7. Two-body loss rate factor as a function of laser power 
for several detunings. 

Figure ^ presents the loss parameter /3 as a function of the 
laser intensity for different detunings 5, from -30 to -5 MHz. 
The uncertainty of the measurements varies from 25 % for 
large detunings to 60 % for small detunings. For all detun- 
ings, P increases with power, which shows that S-P collisions 
are dominant. 
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FIG. 8. Two-body loss rate factor as a function of detuning for 
a fixed intensity / = 80/sat of the laser. 

Fig. ^ shows the loss parameter as a function of detuning 
for a fixed intensity {y~ ~ 80)- For same reason, the rate 
increases when the detuning goes to zero, as the population 
in the P state increases. Our results for /3 agree with previous 
results [p"o|j8| jlll[l^ within the given error bars, extending the 
measurements to a wider range of parameters. For example, 
at -5 MHz and in an intensity range for which /3 is not ex- 
pected to vary strongly (/ = 140 to 200/sat), Kumakura et 
al. [| find P = (4.2±1.2) x 10"* cm^s, Browaeys et al. Q 
fi — 2 X 10"* cm''/s with an uncertainty of a factor 2 and 
Tol et al. /3 = (1.3±0.3) x 10"* cm^/s. Our measurement 
/3 = (3.5±1.4) X 10"* cm^/s agrees best with ||]. One should 
also note that we find neither a decrease of 13 for high intensi- 
ties at small detunings, nor a decrease of P at small detunings 
for a given intensity : this differs from the results of 

0- In 

fact, we find that P increases with intensity at small detun- 
ings, and also increases with decreasing detunings at a given 
intensity. We also disagree with the results of Q where much 
smaller rates are found. 

Finally, we also measured the loss rate in the trapping con- 
ditions {5 — —45 MHz, I = 3107sat) : the decay rate of the 
number of atoms was found to be Pn{0) = 30s"^ at a density 
of 10^'' at/cm^, which gives P = 3 x 10"^ cm^/s. 
One can further analyze these data following the simple model 
of which relates the decay constant P to the constant rate 
coefficients Kss, Kap and Kpp and to the populations of the 
excited and ground state levels, np and tts respectively: 



P = K, 



, TTs + 2Ksp H~ TTp TVp 



(10) 



Experiments or theory have shown that the contri- 

butions KssT^s a-nd KppTTp to the total rate /? are smaller than 
the Ksp term by approximately two orders of magnitude. 
From the measurements of the fluorescence signal in flgure ^, 
we derive np for each experimental point, as F/Fmax in eq. 
(^) is equal to 2 x TTp. 
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FIG. 9. Rate coefficient Ksp for all our measurements, as a 
function of the laser intensity I for several detunings. 

In figure ^ we then plot Kap for the ensemble of our 
data. We do not see clear evidence for a dependence of Ksp 
with the detuning or the intensity within the dispersion of 
our data. To a good approximation, we estimate then that 
Ksp is actually constant in the explored range of parameters: 
Ksp — (1.0±0.4) X 10"^ cm^/s, with a dispersion that roughly 
agrees with the error bars we claim. This result agrees with 
the first measurement ever performed [^, but the precision 
is now much improved. It also agrees well with the measure- 
ments of Q where the authors found Ksp — (8.3±2.5) x 10"** 
cm^/s, assuming that for their parameters {S — —5 MHz and 
7 = 30 mW/cm^), tts = np = 0.5. 

An important point is that, in contrast with the measure- 
ment of the fluorescence at resonance where the transition is 
assumed to be saturated, vTp in the compressed MOT never 
reaches 0.5 in our measurements : even for the smallest de- 
tuning and the highest intensity, np is only 0.2. This explains 
why the results for /3 in figure at 5 = —5 MHz strongly in- 
crease for increasing intensity over the whole explored range. 



VI. CONCLUSION 

We measured the absolute two-body loss rate between 
metastable atoms in a magneto-optical trap as a function of 
detuning and intensity. We extended the range of these pa- 
rameters and compared the results to those of previous mea- 
surements, mostly performed at small detunings. Using a new 
experimental approach, we obtained reliable values for the 
two-body loss rates with an improved accuracy as compared 
to most earlier results. In the region of overlap of parame- 
ters, we find a good agreement with previous measurements, 
within the quoted uncertainties. We find a loss rate monotoni- 
cally increasing as a function of intensity and decreasing with 
detuning. Our measurements are interpreted with a simple 
model, giving a rather constant loss rate Ksp, with an aver- 
age value of (1±0.4) x 10~^ cm'^/s, as already found in the 
very first measurement of jl^]. We believe that the quality 
and the extended range of our measurements should moti- 
vate more theoretical work, in order to understand better the 
peculiar dynamics of Penning collisions between metastable 



helium atoms in the presence of light. 
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APPENDIX A: MODEL OF THE ABSORPTION 

In this appendix we describe the method we used to quan- 
titatively interpret the absorption images of the atomic cloud 
when a standing wave configuration of the probe beam is used, 
and for an arbitrary saturation parameter. We describe the 
atoms as two-level atoms characterized by a non linear sus- 
ceptibility: 



5-i{r/2) 



{r/2)^ + 5^ + \n\y2 



(Al) 



where n{x, y, z) is the atomic density, d the atomic dipole, 5 
the detuning, F the inverse lifetime of the excited state and 
Q. is the Rabi frequency given by 



m 

2 



(+) 



{+), 



+ c.c. 



(A2) 



where E is the electric field. The direction of propagation of 
the beam is z and the field is supposed to be uniform in the 
plane (x,y) . The propagation of the field is then described by 
the Maxwell equations: 



[A + fc^(l + x)]n(^) = 



(A3) 



where feo is the wavevector of the light. 

The principle of the model is to use the slowly varying en- 
velope approximation generalized to the case of a standing 
wave. We then decompose the probe beam field as: 



n{z) =A+{z) e^'^o" -f A_(2)e" 



(A4) 



where A+, A- are the slowly varying envelopes of the wave 
going towards positive z and negative z respectively. A simi- 
lar decomposition holds for the nonlinar susceptibility of the 
atoms: 



x(z) = Xo(^)+X+(^)e'"'"^+X-(^)e" 



+ . 



(A5) 
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where xo, X+ s-id X- are slowly varying envelopes, and where 
we neglect terms in the expansion describing generation of 
frequencies others than the probe frequency via the non linear 
interaction. 



If we insert the expansions ( A4 ) and ( A5 ) into the propagation 
equation (A3) and use the rotating wave appoximation, we 
obtain a set of two coupled differential equations for the slowly 
varying field amplitudes A^. By splitting the complex 
amplitudes into modulus and phase: 



A. 



A^ 



\A\e 



itp- 



(A6) 



and by introducing the phase difference — 0-) in the 
definition of the slowly varying susceptibilities x+ a-nd X-- 



X+ = X+ e 
one can write 



X- = X- e 



d\A+ 
dz 

d\A- 
dz 



= y (Imx+lA_|+Imxo|^+i; 
k 

= -Y(Imx-lA+|+Imxo|A- 



(A7) 

(A8) 
(A9) 



By using expressions (|A1|) and ( [A^ ) , the quantities fco Imx4 
fco Imx- and fco Imxo are readily calculated: 



3X2 

ko Imxo = — n{x, y, z) a fo 



ko Imx4 



3A 

ko Imx- = — n{x, y, z) a fx 



where 



(r/2)2 
1 



; /i = 



\A+\\A-\ 



|A_|2) 
/o 



(r/2)2+<52 + i(|4+|2 + lA_12) 



(AlO) 
(All) 

(A12) 
(A13) 
(A14) 



As a last step we eliminate the atomic density n{x, y, z) from 
the equations by changing variable: 



Ziz) 



n{x, y, z)dz 



and we obtain the final coupled equations: 

d\A, 



dZ 
dZ 



3^' t,^^ 



o\2 

— a{MA+\ + fo\A-\ 



where: 



^- = 4-/(r/2) 



A+ = A+/{r/2) . 



(A15) 

(A16) 
(A17) 

(A18) 



For 5 = and in the limit of small saturation parameters, 
one has q: = 1, /o~l, /i~0 and one recovers the usual 
decoupled equation for l ow s aturat ion a bsorption. We have 
now to solve equations ( A16 ) and (A17). More precisely we 
wish to calculate the column density 



Z°° = 



n{x, y, z')dz' 



(A19) 



for each effective pixel (x,y) of our image of the cloud. For 
each effective pixel, we can measure the initial conditions: 

\A+\^(Z{-oo) = Q,)=h (A20) 
\A- \^{Z{-oo) = 0)=lf (A21) 

corresponding respectively to the intensity of the probe beam 
before passing through the cloud, or equivalently without the 
atoms, and to the intensity of the probe beam that passed 
through the a tomic cloud. For symmetry reasons, the column 
density ( A19 ) is given by 2Zo — Z{0), where Zq verifies 



\A+{Zo)f^\A.{Zo)\' 



(A22) 



For each pixel (x,y), we then integrate equations ( A16 )-( A17) 
numerically using the initial conditions ( A2C )-( p^21[ ) until 
\A+{Z)f = \A-[Z)\^. The corresponding Z multiplied by 2 
gives the column density. Note that, contrarily to what hap- 
pens in the low saturation regime, we here need the values U 
and // separately (and not only their ratio), which implies a 
calibration of our CCD camera. 
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